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. , We have studied the quark-hadron phase transition with RMFT motivated equation of state for a 

strongly interacting hadronic sector and lattice motivated equation of state for weakly interacting 
QQ ' QGP sector. The interactions in hadronic sector are dominated by the exchange of scalar and vector 

mesons (a* a, lj, p, (j>) thereby allowing this phase to be modelled by the interacting baryonic, 
. pionic and Kaonic fields. The pionic and Kaonic fields are incoorporated on equal footing to 

Q-T 

I ■ baryonic field rather than including pions and Kaons as exchange particles only. The effect of 
Oh" 

' interactions on quark-hadron phase transition curve was studied using Gibbs criteria for phase 



cn 
o 

(N 



cquillibrium. It was found that the first order quark hadron phase transition curve ends at a 
critical point, whose coordinates coincide with that of the lattice gauge theory result involving 



> 

f***^ ' lattice reweighting technique. 



I. INTRODUCTION 



Since the discovery of asymptotic freedom \1\ in non-abelian gauge field theories it 
was postulated that a phase transition from nuclear state of matter to quark matter is 



possible. Such transitions were further argued to happen at high temperature and /or 
densities where each participating nuclie would loose its existense to more fundamental 
degrees of freedom i.e quarks and gluons and the new state of matter dubbed as quark 
gluon plasma (QGP) would be formed. Since then a considerable effrort has been put 
forward to create and understand the properties of this new state of matter (QGP) and 
the corresponding phase transition. 

In order to study the phase transition in general one needs to have a complete 
description of a given state of matter on the basis of some underlying theory and only 
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then one can in principle comment on the nature of phase transition by studying the 
variation of some parameter intrinsic to the given state as it goes through the phase 
transition, e.g a study of variation of order parameter in case of Landau theory of 
phase transition [2], which would obtain non-zero value in one phase and vanish in 
other phase, and the manner in which it vanishes would decide on the nature of phase 
transition. 

However such a complete description for hadronic matter on the basis of quantum 
chromodynamics (QCD) is far from being completely understood as it falls in the regime 
where the strong coupling constant is too large for any peturbative expansion to remain 
valid. This opens up an interesting field in strong interaction physics where empha- 
sis is on developing non-perturbative techniques from first principles e.g Lattice QCD 
[3] or studying such non-perturbative systems using some phenomenological methods, 
Hadron resonance gas (HRG) and Bag models [I] etc. However, each of the techniques 
mentioned above have some limited applicability e.g. in case of lattice-QCD one can 
not describe a system with non-zero baryon-chemical potential , therefore it is not 
possible to study the entire QCD phase diagram which involves a correct description 
of hadronic system with large chemical potentials as well ( even though some recent 
advances have made it possible to extrapolate the lattice studies to hadronic systems 
with finite chemical potentials but still systems with large chemical potentials are yet to 
be understood completely). The phenomenological thermodynamical models like HRG 
models and Bag models describe some of the properties of hadron and quark phases 
very well under some simplifying assumptions but by no means offer the complete de- 
scriptions of the these phases. The various lattice and phenomenological studies done 
so far have yielded different results, regarding the nature of phase transition, e.g in case 
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of HRG based model study, one has a first order phase transition from hadronic state 
of matter to quark matter throughout the phase diagram, even though recently it has 
been claimed that after certain point the nature of phase transition changes from first 
to second order [5]. Where as the lattice studies predict the existence of a critical end 
point for this transition, which marks the shift in the nature of phase transition from 
first order to second order or to smooth crossover, depending upon the quark masses 
used in the calculation [6J. Although most lattice QCD calculations indicate the the 
existence of critical end point (CEP) for \xb > 160 Mev [6], however its exact location 
is yet to be established. 

In this work we focus on the problem of quark hadron phase transition, where empha- 
sis is on developing equation of state for strongly interacting hadronic matter very near 
to the phase transition region, using relativistic field theoretical description, in some ap- 
proximaiton scheme, such as a mean field theoretical (MFT) description. It was found 
that for such strongly interacting hadronic matter the interactions are dominated by 
the exchange of certain mesons (a*, a, u, p, 0). Therefore one could in principle model 
the entire hadronic phase in terms of interacting baryonic, pionic and Kaonic fields, 
where pions and Kaons are treated on equal footing to the baryonic field rather than 
considering these pions and Kaons as exchange particles. This scenario will supposedly 
remain more and more valid as we approach more and more closer to the phase tran- 
sition region, which happens to be our main region of interest for studying the quark 
hadron-phase transition. Quite clearly as we move away from such regions our simpli- 
fying assumption will breakdown and 7r meson and K meson exchange among various 
baryons has to be taken into account. For QGP Sector we use an equation of state that 
takes into acount the perturbative interactions among quarks and is consistent with 
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lattice data. 

The remainder of this article is organized as follows: In section II we present the 
formalism used in this work. In section III we show the numerical results and discussion. 
Finally in section IV we summarize the results and give the brief concluding discusion. 



II. THE FORMALISM 

We present equation of state (EOS) for the hadronic phase and the QGP phase used 
in this work along with the definitions 

A. Quark-gluon phase: Quarks u,d and s (+ gluons) 

To study the QGP phase with three quark flavours (u,d,s) and gluons, we use a Bag 
model [1] equation of state, with perturbative corrections of the order of a s [7] . The 
pressure, energy density take the form [8] 

x [/,,/(*) +7,,/(*)] -b (i) 

/ A \ N r°° i Nf N t°° i 

e(T,n f ) = (l--t)^Lj o k 2 dk^ + ml f g (k) + (1 - C) J2 ~$ J k 2 dk^V+ m f 

x [/,,/(*) + /,,/(*)] +B (2) 

where ( = a s is a model parameter and 1 — £ , 1 — |£ represent the perturbative cor- 
rections to the kinetic terms of quarks and gluons, respectively. The non-perturbative 
vaccum effects are taken into account by the using a Bag constant B. The m g , m q 
denote the gluon and quark masses and f g (k) , f q (k) are the bosonic and fermionic 
distribution functions, respectively. Here N g = 2(N% — 1) is the number of transverse 
gluons and N c is the number of quark colors (iV c = 3). In case of two massless quarks 
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(u,d) and a massive (s) quark, alongwith gluons, the above expression for the pressure 
and energy density reduce to: 
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x(E 2 ~m 2 s )i(f k + f k )+B 



(4) 



here iV 9 = 16(1 — ~£) is the "effective" number of the gluons and Nf = 2(1 — £) is the 
"effective" number of light flavours. However it should be mentioned that throughout 
our discussion we use QGP phase with zero net strangeness. This assumption seems to 
be very reasonable in case of heavy ion collision, where the colliding nuclei do not carry 
any strangeness and therfore no net-strangeness is carried over to the QGP phase and 
the new thermally generated strange quark- antiquark pairs in the QGP phase always 
maintain the zero net strangeness condition. 

B. Hadronic phase: baryons, Kaons and pions 

In nuclear physics where there are several types of quanta of the nuclear field, dif- 
fering in quantum numbers and masses (in contrast to the electrodynamics with only 
one field quanta), the nucleon-nucleon (NN) potential is defined as the superposition 
of components with different space time transformation properties and different radii 
of action. In the calculation of the NN potential the range of potential is divided into 
three parts, the external part r > 2.14/m, the intermediate part 0.71/m < r < 2.14/m 
and the internal part r < 0.71/m [9]. These three regions are determined by different 
mesons and resonances that make up the contributions to the NN potential (from the 
uncertianity principle it follows that the radius of the action of the nuclear force is 
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determined by the compton wavelength of the respective quanta). At present it can be 
considered as an established fact that the external and intermediate parts of the NN 
potential in meson theory are determined by the meson exchanges, however the internal 
part, which corresponds to the internucleon distance less than 0.71fm, this picture of 
the single boson exchange is not acceptable, as many processes make contributions to 
the NN forces in this region. Therefore this part of the interaction is mainly described 
by the phenomenological methods or by employing quark-hadron dynamics. The ex- 
ternal region is completely determined by the one pion exchange. The pion may be 
coupled to the nucleon field by either pseudovector or pseudoscalar coupling. How- 
ever sometimes a mixture of both coupings is used [9]. To describe the intermediate 
range of NN potential a proper treatment of two pion exchange contribution is required 
along with the exchange of other mesons. However in many NN models the attrac- 
tion produced by 2ir exchange is simulated by scalar-isoscalar a mesons with the mass 
around 500-600 MeV [9]. Keeping in view the fact that for a possible hadron-quark 
phase transition, hadrons must come close to each other sufncent enough so that the 
quarks from each individual hadron enter the state where quarks are asymptotically 
free, therefore the inter-nucleon distance cannot correspond to the external part of the 
NN potential (r > 2.14/m). Thus the one pion exchange which describes the exte- 
rior part of NN potential can be safely ignored in comparison to the scalar and vector 
mesons (a*, a, u, p, <fi) exchanges. However since the production of pions and Kaons 
is large at finite temperature (T), hence we need to incorporate these in the system 
along with the other hadrons (not as exchange particle though) with their intercations 
to baryons governed by the exchange of scalar and vector mesons. 

Therefore in order to describe the hadronic phase we develop a relativisic field the- 
oretical model in which baryons, pions and Kaons are included on same footing, and 
are interacting via the exchange of mesons. The baryons considered are (N, A, S, S, A) 
and the exchange mesons include isoscalar-scalar and vector mesons (er, u), isovector- 
vector mesons (p) and two additional hidden strangeness mesons (a*, 0) The lagrangian 
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for the hadronic phase can therefore be written as, 



c Total =£ B+ Cr + £n ( 5 ) 

Where the effective lagrangian for baryons is [10] : 

Cb = ^ ^ B _ m B - 9<jBO + 9a' BO* - 9ujB1^ ~ g^Bj^ - 9 'pBl '^Tj ]* B 

B 

+ \d^a - l -mla 2 - l -g 2 a* - \g z a A - \w, v W^ + ^m^w" 

- \S» V S» V + \m%^ (6) 



where = d^ u - d"^, R afny = d»p av - d" p a » + g p t abc p b ^p cu and W^ v = d^uj v - d u u» with 
the sum carried over the entire baryon octet. For Kaons [TT] the effective lagrangian 
takes the form: 

C K =J2 d v I<] d ^ K - m t K] K - g*Km K K ] Ka - g a » K m k K ] ' Ka* - g uk Kk1%Ku li 

K 

- gpkK^TKU^Kp' 1 ~ g^Kkd^K^ + (g^KUp + g P Krpp + g^^fK^ K (7) 
Now similarly for pions one can write 

Ar = d^d^n - mlnkr - g^m^Txkra - g un %kd p Ttu) fJ ' - gp^T^i^np'" 1 

7T 

+ (Pwtt^ju + gp^Ppf IT^IY (8) 

One can easily see that both of these coupling schemes mentioned above fulfill the 
Ward identity. The resulting field equations for fermions and bosons obtained after 
minimising the corresponding action, S = J Cd 4 x, are non-linear in form and because 
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of strong coupling constants cannot be solved using perturbative techniques and hence 
are to be solved using an approximation scheme e.g in mean field theoretical description. 
Using such an approximation, the field equations for the mesons (a, u>, p, a*, 0) turn 
out to be as follows: 
For sigma (a) field; 
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d 3 k 



-.m* B (n B + n B ) 



(2tt) 3 J 2uj b 



(n b + n b ) 



(9) 



similarly for omega and rho fields one has, 



m 2 „uj - 



C3^ 3 =^9uB 7^5 J d 3 k(n B -n B ) + 2g ub -^-^ J ^ (E£ n b + E b n b ) 



b=(K,Tr) 



b =(K,-K) 



Ei 2 \ 7 f 

(2^.9^ + Ig^bgpbpTz + 2g ujb g^) x -^—^ J — (n b 



■ + n b ) 



(10) 



m 2 p P = Y,9 P BT3 j d 3 k(n B -n B ) + ]T 2g pb T 3 J ^ (E+ n b + E, 



b =(K,ir) 



E ( 2 9^ b 9pbUT 3 + 2g 2 pb rip + 2gt >b g pb T 3 $) 

b=(K,ir) 



7 f d 3 k 



(2tt) 3 J 2iv b 



(n b + n b ) 



(11) 



For sigmastar (a*) field we obtain, 
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b=(K,TT) 



El f d 3 k . _ . 

(12) 



and finally for the phi field we have, 
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J d 3 k(n B - n B ) + 2g$ h -^-^ J ^ (E b + n b + E b n b ) 



b=(K,w) 



^2 { 2 9^bg<i>bU + 2g ib g pb T 3 p + 2$g| b ) x -^-^ J ^ (n b + n b ) 



b=(A» 



(13) 



with the distribution function for baryons and antibaryons given by: 

n B = [exp(E B -v B /KT) + l}- 1 
n B = [exp (E* B + Vb/KT) + If 1 

where the effective mass and chemical potential for baryons is, 



(14) 
(15) 



m* B =m B +g aB a + g<j* B a* 

v B * =p B - g^ B Lo - g^ B 4> - Q p bTzbP 

for Kaons and pions the effective mass takes the following form, 



(16) 
(17) 



m* k =^J m\ + m k (g ak o + g^-kcr*) 



(18) 
(19) 



and the effective chemical potential: 



vk =Mfc - g^ku - g<j>k<i> - gpkT-ip 
=/x w - g.^to - g p7T r 3 p 



(20) 
(21) 



The pressure for the strongly interacting hadronic matter can be derived using the 
energy momentum tensor given by 
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rptotal ptotal , ff) (f\ \ 



total 



Q£total 



(22) 



as 



ptotal _ . rptotal . 



(23) 



Where indices (i) on is summed over and each component corresonds to different 
field. $*=(\I/,K,7r ), where (using Eq.5) 



rptotal _ rpbaryons , rpkaons , rppions 
ii ii ' H ' ™ 



(24) 



Thereby allowing pressure to be written as a sum of three parts 



ptotal pbaryons pkaoris _|_ ppions 



(25) 



where the pressure due to baryons and kaons turns out to be, 



k dh . . 1 o o ^ 3 ^ 4 ^ 2 2 ^ 4 

(n fc + ra fe ) - -m a a - -g 2 a - -g 3 a + -m^uj + -c 3 uj 



id , 



B 

,1 22 1 2 *2 . 1 2/2 



(26) 



jKaons 



1 v-^ 7 /" " & > 2 / - x 1 W/ , \2 * 

3 2^(^)3 y ^~ fc {n k +n k ) + -2_^{{g UJK uj + g (t>K (j> + g pK Tz P ) -g a > K m K a 

\ y^( 2 9uKU + 2g ujK uj + 2g pK r 3 p 



gaKm K cr - m 2 K ) 



2g<t>K4>) 



(2 



7 f d 3 k . 
WJ^ {nk + nk) 



7 f d a k , , 
?r) 3 7 2w fc 



K 



(27) 
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Similarly the contribution to the total pressure from pions comes out to be, 



ppions = 1 J- J ^A k * („ fc + n k ) + 1 + g p7T T 3 pf - g^m^a - ml) 



d 3 k 



1 v-^ r 7 f d 3 k 

- -22( 2 9^ + 2 9p^T3P) j—^; j —(E+n k +E-n k ) 



(28) 



Now the energy density for this strongly interacting hadronic matter can be calculated 
as follows, 



^.total ^^>| . rptotal . 



total 



(29) 



which after using Eq.5, simplifies to, 



^total ^baryons _|_ ppions ^.kaons 



(30) 



where the energy density for baryons is found out to be, 

£b = ^2 J~y J fc2 \A 2 + m *B 2 ( nk + n k )dk + ^m a 2 a 2 + ^92<? 3 + j93^ 4 + ^ 



2 2,° 4 
77J W OJ + ^C 3 W 



i ^ 2 2,1 2 *2 . 1 2/2 



(31) 



similarly for Kaons and pions one has, 



7 f d 3 k 
2u k 



£K =Hj<^Y J 2^£ [ (jB+)2 Uk + ^ ) 2 M+( m k 2 +9vKm K a + 9 a *Km K a*)j^ J 

, _ , , x 2 7 f d 3 k _ 

x (n fe + n k ) - (^ifw + g^K + g P KT 3 p) j—^ J — (n fe + n fe ) 



(32) 
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2 \ 7 f d 3 k 



(33) 



In this work we will use a successful parameter set of RMFT model, TMI (Table-I) 
these parameters have been determined by fitting to some ground state properties of 
finite nuclie, including unstable nuclie. As for meson-hyperon couplings we take the 
naive quark model values for the vector couplings following [10] (Table-V). However 
regarding the antibaryon couplings, there is no reliable information suitable for the 
high density matter and therefore antibaryon-meson coupling constants motivated by 
G-parity transformation will be used (Table VI). It is important to note that this simple 
consideration based on the G-parity transformation of meson fields is certainly an ide- 
alization as there are several effects in many body systems that can distort this picture 
and one is forced to use modified antibaryon-meson couplings [12] . On phenomeno- 
logical level this can be taken into account by multiplying the baryon-meson couplings 
with a modifying factor rj which will act as a free parameter in the model and whose 
value can be varied from < rj < 1 to take control of maximally strong antibaryon 
couplings to non-interacting antibaryons. However this requires some detailed study 
and we will not take this into account here. The coupling constants for Kaons and 
pions will be as in Ref. [32] and this completes our description for the hadronic phase. 

III. RESULTS 

First of all we will present some of the features of the QGP phase followed by the 
description of hadronic phase and associated phase transition along with other observ- 
ables. In case of QGP phase, the EOS contains two undetermined parameter's B and £ 
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which are to be fixed uniquely for making any valuable prediction. As far as bag value 
B is concerned, we fix the bag value to be B = 344MeV/ fm 3 following Ref. [H], |15j . 
where it was successfully used for the fluid-dynamical calculations of the heavy-ion col- 
lisions. However in view of the fact that bag value can itself be a function of chemical 
potential and temperature, keeping the bag value B fixed, is an approximation. We 
will stick with this approximation as the more general form of B = B(fi,T) requires 
a detailed study and will be picked up in future studies. Now to determine the model 
parameter ( we first see its influence on the QGP equation of state. Fig(l) shows the 
variation of pressure with temperature for different values of paremeter (. Here ( = 
corresponds to the non-interacting QGP phase. To choose among the possible values of 
the parameter (, which can vary in the range [0-0.3] [16] . we show in fig(2) the scaled 
energy density of QGP phase for B = 344MeV/ fm 3 along with the lattice calculation 
results. It can be easily seen that the the best fit occurs for the ( = 0.2 




-400 1 1 1 1 1 1 1 

50 100 150 200 250 300 

T (MeV) 

Figure 1: Pressure in the QGP phase as function of temperature for £ = 0.0, 0.1, 0.2, 0.3 for a, b, c and d 
respectively. Here fi q =0 MeV. 
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Figure 2: Scaled energy density for QGP phase with £ = 0.0, 0.1, 0.2, 0.3 for a, b,c, d respectively and To is 
the critical temperature corresponding to each value of £ [8]. Solid circles are for lattice data set with Nf = 3 



Throughout our calculation we will use these values of bag constant (B) and model 
parameter Now with these fixed values of bag value (B) and model parameter ( we 
next plot in Fig(3) the variation of pressure with temperature (T) for different values 
of chemical potential \i q . 




50 100 150 
T (Mev) 



200 



Figure 3: Pressure as function of temperature for different values of chemical potential, B = 344MeV//m3, 
C = 0.2, fi q = 220, 190, 170, 150, 135, 120 MeV for a, b, c, d, e and f, respectively 



Next we show some of the features of the hadronic phase. The parameter set we use 



3 In present model the parameter f needs to be determined by comparing scaled enery density and lattice data at 
finite fig , however in abscence of lattice data at finite fig this value of f will serve as an approximate value 
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here are listed in table (I, II, III and IV). In Fig(4) we plot the variation of pressure in 
this phase. 
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Figure 4: Pressure for hadronic phase with /i 9 = 120, 135, 150, 170, 190, 220 MeV for a, b, c, d, e and f, respec- 
tively 

The Pff-T curves develop the cusp which becomes more and more evident as one de- 
creases the chemical potential \i q . This behaviour can be attributed to the interplay 
between attractive and repulsive attractions present in the system. Now if we calculate 
the slope of any general curve with this feature as in Fig. 5a we conclude that, for the 
region AB the entropy desnity is: 



large therefore the probability of finding the system with the pair of values (fj, q , T), corresponding to 
region AB is given by: 




(34) 




(35) 
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where Sab = 



—S and ks is the Boltzman's constant. For the point B we have: 



sab =0 



(36) 



which imply 



s = k B \nn = 



(37) 



as a result of which £1 = 1, where are the number of states accesible to the system. 
Also = E- f , where f are the number of degrees of freedom for the system under 
consideration. Therefore the point B corresponds to the configuration where the number 
of degrees of freedom are very small infact zero. This configuraion can be therefore 
thought to be that of a system, wherein particles under the influence of strong attractive 
and repulsive forces attain the equilibrium state and arrange themselves in lattice like 
network. It is clear from the arguments above that the probability of finding the 
hadronic system with chemical potential and temperature corresponding to the region 
AB is very small as compared to the region BC and beyond. Therefore the effective 
P H — T curve can be drawn as in Fig. 5b. 
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Figure 5: Py-T curve for hadronic phase for /U g =0 MeV 
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Now to describe the quark-hadron phase transition we apply Gibbs criteria of phase 
transition. This hypothesis requires following set of relations to be valid at the phase 
coexitence points [T7] . : 



Tq=T h 

where it is required that the state with higher pressure is more stable. Now, following 
three cases can arise as are shown in Fig.6. Here we show the intersection points be- 
tween the P-T curves of hadronic and QGP phase. Out of the these intersection points 
only first two, correspond to the phase coexistence among the QGP and hadronic phase, 
however for intersection point in case 3, the pressure due to quarks and gluons is still 
less than bag pressure B, therefore the QGP phase corresponding to this intersection 
point is not stable and therefore cannot coexist with the hadronic phase. 




60 80 100 120 140 160 180 200 220 60 80 100 120 140 160 180 200 220 60 80 100 120 140 160 180 200 220 
T(MeV) T(MeV) T(MeV) 



(a) fi q = 135 MeV (b) fi q = 120 MeV (c) ^ q = 100 MeV 



Figure 6: P-T curves for hadronic and QGP phases. For lower values of chemical potential the intersection 
point starts to appear below the P=0 axis, where the QGP phase is unstable. 
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Now using the same procedure for the entire set of baryon-chemical potential \ib values, 
the quark-hadron phase diagram turns out to be as in Fig. 7 
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Figure 7: quark-hadron phase transition curve 

The first order phase transition curve ends up at the critical end point CEP (ci) (/x 9 =120 
MeV, T=154 MeV), which is close to the critical end point point as found using lattice 
Resummation technique (LR2) [18]. In fig.7 we also show the critical end points as found 
in different studies [19]. Next we plot the variation of net baryon density (tib) along this 
curve. Fig.8 shows the variation of ng with critical baryon chemical potential Hb c - For 
the sake of comparison we also plot the the variation of net baryon density % in case 
of hadron resonance gas (HRG) model. It is interesting to see that the number density 
in case of present model saturates beyond a particular chemical potential, and this 
saturation happens in between the HRG curves with r=0.8fm and r=0.7fm. Therefore 
it can be assumed that the time independent scalar and vector potentials force the 
hadrons to attain an effective size. This makes sure that the number density of hadronic 
phase with point particles does not rise high enough to make hadronic phase stable at 
very high temperature or chemical potential. Next we show the variation of % with 
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Critical temperature T in Fig. 9, here the solid dot corresponds to CEP (ci). 
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Figure 8: net baryon density v/s critical baryon chemical potential, for RMFT and HRG (with excluded 
volume assumption) modelled hadronic phase, where r=0.6 fm,0.7fm and 0.8 fm for c,d,b respectively 
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Figure 9: mall net baryon density v/s critical temperature, for RMFT and HRG, here r=0.6 fm, 0.7 
fm, 0.8fm for c,d and b respectively. 



To elaborate this point further, we plot the variation of proton density with temperature 
in our present model and compare it with HRG model for three different hard core 
volumes, this is shown in Fig. 10. 
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Figure 10: Proton density (np) v/s Temperature (T) for chemical potential fj.B = 405 Mev, r=0 fm, 0.683 
fm, 0.8 fm for a,d and c respectively. 

The proton density does not rise to an arbitrary large value in our model, but saturates 
to a value of 0.06/m~ 3 , this is exactly the saturation as found in HRG calculation with 
hard core radius r=0.683fm. Next we calculate the response of the hadronic medium to 
the external perturbations, for this we calculate the velocity of sound in the medium, 
whoose square at vanishing chemical potential is given by [2D] 



In Fig. 11 we plot the variation of c 2 with energy density, followed by Fig. 12 where we 
show variation of c 2 with temperature. 



2 dP 



(44) 



since, P = P(T,/i = 0) the above equation can be written as 




(45) 
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Figure 11: <? 3 -e for a hadronic phase with a chemical potential jj, q = MeV 
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Figure 12: c 2 g -T for the hadronic phase with a chemical potential fi q —0 MeV 

Starting with a value close to the ideal gas limit cj. — > l/3,there appears a sharp dip 
where c 2 goes to a minimum value of 0.008, beyound which it starts to rise again but not 
sharply though. This softening is usually attributed to the possible phase transition, 
however keeping into account the dynamics of present model, where scalar and vector 
interactions force particles to attain a state with degrees of freedom T equal to zero, 
this temperature 'T' for which c| dips corresponds to this very configuration where 
we expect, given a perturbation the disturbance generated propogates through small 
extension of hadronic system rather being transmitted and hence a value c 2 = 0. It is 
only after this state a possible phase transition can occur. 

IV. SUMMARY 
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We have presented a study of deconfinement phase transition from hadronic matter to 
the quark gluon plasma (QGP) that could be formed in the heavy ion collision. We 
modelled the entire hadronic phase with pions,kaons and baryons, with interactions 
among various hadrons carried by some specified vector and scalar mesons only. The 
QGP phase with perturbative intercations, was modelled in terms of MIT bag model 
with an equation of state consistent with the lattice data. In such a hadronic phase 
we found that the interactions among various hadrons force each hadron to pick up an 
effective volume, the corresponding radius is found out to be r=0.683fm. This has a 
serious consequence as it makes sure that the number density of hadrons does not rise 
to an arbitray large values, which makes the hadronic phase stable again at very large 
temperature or chemical potential. We further found that a first order phase transition 
is not possible throughout the phase diagram and ends up at a critical end point (CEP) 
c\. This is found to lie very close to the CEP as fouund in LR2[18]. However it should 
be mentioned that at higher temperature and/or chemical potential one needs to take 
into account the fluctuations about the mean value also into account, which calls for 
an approach beyound the mean field level, this surely needs a systematic study and is 
currently under investigation. 
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Table I: The parameter set TMI [10] used in the calculation. The masses are given in MeV 



m N 


m a 


m w 


m p 


9aN 


guN 


9pN 


92(fm x ) 


#3 


c 3 


938.0 


511.198 


783.0 


77.0 


10.029 


12.614 


4.632 


-7.2323 


0.618 


71.308 



Table II: Hyperon coupling constants 





5<x£ 




9<y'A 


g<7*~ 


6.170 


4.472 


3.202 


7.018 


12.600 



Table III: Kaon coupling constants: 



9<jK 




9pK 


9a* K 


9<t>K 


1.93 


3.02 


3.02 


2.65 


4.27 



Table IV Pion coupling constants: 



9np 




97TO- 


9^(7* 


9tt<p 


0.506 


-0.001 


-0.170 


0.0 


0.0 



Table V: Vector Couplings used in the present calculations [TO] 

~^9loN 2~9u}h ~^9u<j 9u>S 
9pN = \g P cr = 9pS, 9pA = 0. 
2<?0A = 2g<f>a = 2g ( j,s= -^j^-guN, g<j)N = 0. 

Table VI: G-parity motivated antibaryon coupling constants [12] 

g^B = g°B, g U B = -gu>B, g P B = gpB- 
2g<pA = 2g<i>Y, = g<ps = ^-^f^g^N g^u = ~9ojn 
9<fB — -g<t>B 
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